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1 .  INTRODUCTION 


Lightweight  composite  materials  have  been  used  for  many  years  in  the 
aerospace  industries.  Some  of  the  most  needed  design  informations  are  the 
failure  mechanism  and  the  failure  criteria.  Determination  of  the  failure 
mechanism  and  the  failure  criteria  requires  in  general  a  rigorous  stress 
analysis  at  the  singular  point  in  the  material.  The  lack  of  such  an 
analysis  is  one  of  the  impediments  in  understanding  the  failure  mechanism 
and  in  determining  the  failure  criteria.  In  the  case  of  layered  composite 
materials  in  which  each  layer  is  of  orthotropic  material,  experimental 
observations  indicate  that  the  failure  modes  are  in  general  either  along 
the  interface  between  the  layers  or  transverse  to  the  layers.  For  instance, 
consider  the  layered  composite  shown  in  Fig.  1  in  which  each  layer  is  a 
fiber  reinforced  composite  laminate.  The  fibers  lie  in  the  plane  of  the 
layers  although  the  orientation  9  of  the  fibers  may  vary  from  layer  to 
layer.  When  the  composite  is  subjected  to  an  extensional  strain  in  the 
Xj-direction,  a  delaraination  may  occur  along  the  free-edge  MN.  The  stress 
singularities  at  the  free-edge  such  as  the  point  M  were  investigated  in 
[1].  It  was  shown  in  [1]  that,  except  for  certain  special  combinations  of 
the  ply-angles  on  both  sides  of  the  interface,  the  stresses  have  the 
logarithmic  singularity  at  the  free-edge  in  addition  to  the  r  (k  >  0) 
singularity.  Moreover,  unlike  the  r  singularity  whose  existence  depends 

on  the  stacking  sequence  of  the  layers  and  the  complete  boundary  conditions, 
the  existence  of  the  logarithmic  singularity  at  point  M  depends  only  on 
the  ply-angles  on  both  sides  of  the  interface.  Instead  of  the  delamination, 
a  transverse  crack  shown  by  the  dashed  line  may  occur  in  the  layer.  There 
are  stress  singularities  along  the  transverse  crack  edge  MQ  which  is  on 
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the  interface  between  two  layers  and  hence  may  initiate  a  delamination  along 
the  interface.  The  nature  of  the  singularities  at  the  crack  edge  such  as 
point  0  is  the  main  interest  of  this  investigation. 

For  this  purpose,  we  take  point  0  as  the  origin  of  the  rectangular 
coordinate  system  in  which  the  x^  =  0  plane  is  the  interface 

while  the  transverse  crack  lies  on  the  x^  =  0  plane.  Fig.  1.  In  the 
(x^,x_)  plane.  Fig.  2,  the  x_-axis  is  the  interface  between  the  two  materials 
and  a  crack  in  material  1  is  located  along  the  negative  X2-axis.  The 
deformation  at  and  near  the  origin  0  is  certainly  three-dimensional,  i.e., 

the  displacements  u^  are  functions  of  x^,  x,  as  well  as  x^ .  However, 
it  was  shown  in  [2]  that,  to  the  first  order  of  approximation,  the  order  of 
singularity  at  point  0  is  affected  by  the  dependence  on  x^^  only  to  the 
extent  of  the  strain  component  at  point  0.  Since  the  problem  is 

linear,  we  will  first  study  the  singularity  at  point  0  due  to  the  dis¬ 
placement  field  in  which  u^  are  functions  of  x^  and  x^  only.  We  will 
then  study  the  singularity  at  point  0  due  to  a  uniform  extension  in  the 
Xj^-direction  and  see  if  additional  singularities  are  present. 

IVhen  the  materials  on  both  sides  of  the  interface  are  isotropic,  the 
singularity  at  the  tip  of  a  crack  which  is  normal  to  and  ends  at  the  inter¬ 
face  has  been  studied  by  many  investigators  (see  [3-8],  for  example).  In 
real  composites,  each  layer  is  a  fiber  reinforced  laminated  material  and 
hence  should  be  regarded  as  an  anisotropic  material.  The  problem  of  a  crack 
terminating  at  the  interface  of  two  anisotropic  materials  does  not  seem  to 
have  been  investigated,  although  the  case  of  orthotropic  materials  was 
studied  in  [9,10].  In  Section  2  we  present  the  basic  formulations  for  aniso¬ 
tropic  materials  in  two-dimensional  deformation.  Instead  of  using 
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Lekhnitskii's  approach  [11]  which  breaks  down  in  degenerate  cases  [12]  and 
requires  a  special  treatment  for  the  degenerate  cases,  we  employ  the 
analysis  which  was  originally  due  to  Stroh  [13]  and  further  developed  by 
others  [14,15]  for  studying  surface  waves  in  anisotropic  elastic  solids. 

In  Section  3  we  study  the  stress  singularities  at  a  crack  tip  in  a  homo¬ 
geneous  anisotropic  elastic  material,  i.e.,  the  material  is  not  layered. 

The  same  problem  was  studied  in  [16,17]  for  anisotropic  elastic  solids 
whose  material  property  is  symmetric  with  respect  to  the  =  0  plane,  and 
in  [is, 19]  for  general  anisotropic  materials.  By  a  different  approach 
from  that  of  [16-19],  we  show  that  K  =  0.5  is  a  triple  root  singularity 
for  a  crack  tip  in  a  homogeneous  anisotropic  medium.  In  Section  4  we 
formulate  the  equations  for  determining  the  order  of  singularities  at  the 
tip  of  a  crack  which  is  normal  to  and  ends  at  an  interface  between  two 
general  anisotropic  materials.  We  also  show  how  one  can  determine  the 
number  of  roots  <  in  a  given  region  in  the  complex  K-plane.  Numerical 
examples  are  presented  in  Section  5  in  which  both  materials  across  the 
interface  are  of  the  same  orthotropic  materials  although  the  orientations 
of  the  axes  of  symmetry  are  different.  It  is  shown  that  there  are  three 
positive  roots  for  k  and  that  all  three  have  the  values  between  0.333 
and  0.698.  The  fact  that  we  have  three  positive  k  and  that  they  are  near 
the  value  0.5  is  not  surprising  in  view  of  the  results  obtained  in  Section  3 
for  the  crack  in  a  homogeneous  anisotropic  medium.  Finally,  in  Section  6  we 
study  the  stress  distribution  near  the  tip  of  the  crack  which  is  normal  to 
the  interface  subject  to  a  uniform  extensional  strain  in  the  x^^-direction. 

We  show  that  a  uniform  stress  solution  exists  and  hence  no  additional  stress 
singularity  is  present  due  to  the  uniform  extension. 
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ANISOTROPIC  MATERIAL  IN  'nVO-DIMENSIONAL  DEFORMATION 

In  the  rectangular  coordinate  system  (x^,x,3,  let  the  displacements 

u.  (i  =1,2,3)  and  hence  the  strains  £. .  and  the  stresses  a.  .  be  functions 
1  13  ij 

of  x^  and  x^  only.  The  strain-displacement,  stress-strain  and  equilibrium 
equations  can  be  written  as 

.  =  (3u./3x.  +  3u./3x.)/2  (1) 

13  13  3 


a.  .  =  c.  £, 

13  i3kp  kp 


3a._/3x_  +  3a., /3x,  =  0 
’  ^  ^  10  3 


where  repeated  indices  imply  summation  and  . 

‘^ijkp  ■■  ^jikp  ~  ^kpij 

are  the  elasticity  constants.  A  general  solution  for  Eqs.  Cl. 2, 3)  can  be 
obtained  by  letting 

u.  =  u.f(Z)  C 

11 

Z  =  x^  +  px  I 

where  p  and  are  constants  to  be  determined  and  f  is  an  arbitrary 

function  of  Z  [13-15].  Substituting  into  Eqs.  (1,2,3)  we  obtain 

a. .  =  T. .  df/dZ  ( 

13  13 


where 


“ik“k  ■  “ 


^ij  '■^ijkZ  ^  ^^ijk3^^k 


""ik  =  "iZkZ^P^'^iZkZ^^iZkZ^  '^iZkZ 


For  a  non- trivial  solution  of  u^,  it  follows  from  Eq.  (8)  that  the  deter¬ 
minant  of  D.,  must  vanish.  That  is, 
ik 
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(11) 


This  results  in  a  sextic  equation  for  p.  Since  the  eigenvalues  p  are  all 

non-real  [11,13],  there  are  three  pairs  of  complex  conjugates  for  p  which 

will  be  denoted  by  p,  and  p  (L=  1,2,3)  where  an  over  bar  denotes  the 

complex  conjugate,  and  three  pairs  of  associated  eigenvectors  u.  and 

1 ,  L 

u.  (L= 1,2,3).  It  should  be  pointed  out  that  some  elements  of  t  given 
1 9  L  1  j 

by  Eq.  [9)  ^^e  related  to  each  other.  We  can  show  from  Eqs.  (4),  [8)  and 
flO)  that  T. .  =  T..  and 

^  ij  ji 

"12  =  -P^13  , 

(12) 

^22  =  -P^23  =  P''^33 

They  agree  with  Lekhnitskii ' s  results  in  [11]  except  that  is  taken  as 

unity  in  [11]:  However,  x^^  vanish  in  degenerate  cases  [12]  and  hence 
Lekhnitskii 's  approach  requires  a  special  treatment  whenever  x^^  “  9* 

The  general  solution  for  the  displacements  and  stresses  can  now  be 


written  as 


“1  ’  ^ 


where  the  summation  is  over  L= 1,2,3  and  f  and  g  are,  respectively, 

L 

arbitrary  functions  of  and  Z^.  Notice  that  Eq.  (13)  contains  six 

arbitrary  functions.  IVhen  p  has  a  double  root,  the  solution  given  by 
Eq.  (13)  is  not  general  unless  one  has  two  independent  associated  with 

the  douple  root  p.  Similarly,  Eq.  (13)  may  not  be  a  general  solution  when 
p  is  a  triple  root.  A  modified  solution  which  maintains  its  generality 
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when  p  has  a  multiple  root  was  given  in  [20].  For  is'  tropic  materials, 
p  =  i  and  is  a  triple  root.  However,  one  of  the  triple  roots  is  associated 
with  the  anti-plane  deformation  while  the  other  two  which  can  be  considered 
a  double  root  are  associated  with  plane  strain  deformations  [21  ] . 

To  study  the  singularities  at  a  crack  tip,  we  let 


‘  vi'7 


(L  not  summed) 


(15) 


where  <,  A  and  B  are  arbitrary  constants  which  may  be  complex.  Intro- 
L  Ij 

ducing  the  polar  coordinates  (r,<J)): 


=  r  cos<})  , 


X,  =  r  sin4) 


(16) 


Equations  (13)  and  (14)  can  be  written  as 


where 


a.  .  =  r 

L 


=  coscj)  +  Pj^sinc}) 


?{Vij,L^L  ^®L'^ij,L^L  } 


(17) 

(18) 

(19) 


We  see  that  if  the  real  part  of  <  is  positive,  is  singular  at  r  =  0, 

When  <  is  real,  one  may  choose  without  loss  of  generality  [1] 


(2Q) 


where  a  and  a  are  real  constants.  Equations  (17)  and  (13)  then  have 

L  Li 

the  real  e.xpressions 


(21) 
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ij  ^  \  L  V  ij.L^L  /  L  \  ij,L^L  // 

wliere  Re  and  Im  stand  for  real  and  imaginary,  respectively. 


where 


may 

also  use 

c.  . 
ij 

instead 

of  c . 
1 

and 

jkp 

a. 

=  C.  .  , 

c .  .  = 

c .  . 

1 

iJ  3 

Jl 

a 

=  0 

a. 

=  0  „  . 

0, 

=  a 

1 

11  ’ 

L, 

22 

33 

a 

=  0  , 

=  0  ,  , 

a 

=  0  „ 

4 

23  ’ 

5 

13  ’ 

6 

12 

£ 

=  £ 

=  £  »  . 

£ 

=  £  - 

1 

11  ’ 

2 

22  ’ 

3 

33 

"4 

= 

“  23  ’ 

■^13  ’ 

^6 

^^12 

We  will  also  write  the  inverse  of  Eq.  (23)  as 


e.  =  s. .a.  , 
J-  iJ  3 


s. .  =  s  . . 
^3  Ji 


where  s^^  are  the  elastic  compliances, 


(24a) 


(24b) 


5.  SINGULARITIES  AT  A  CRACK  TIP  IN  A  HOMOGENEOUS  ANISOTROPIC  MATERIAL 

In  this  section  we  assume  that  material  1  and  material  2  in  Fig.  2 
are  the  same  material.  Therefore,  there  is  no  interface  and  the  crack  is 
in  a  homogeneous  medium.  This  problem  has  been  studied  in  [16,17]  where 
the  material  property  is  assumed  to  be  symmetric  with  respect  to  the 

=  0  plane.  No  symmetry  is  assumed  in  this  section.  Sih  and  Chen  [iS] 
and  Hoenig  [l9]  also  studied  the  problem  without  assuming  the  symmetry, 
but  the  analysis  presented  here  will  show  that  <  =  4  is  a  root  of  multi¬ 
plicity  three.  The  boundary  conditions  at  the  crack  surface  are 

=  0  at  1)  =  tTT  (26) 
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From  Eq.  [19)  we  have 


±7Ti 

e 


at  (J)  =  ±iT 


(27) 


and  application  of  Eq.  (26)  to  (18)  yields 


-K7T1  „  -  Kiri) 

e  +  B,  T, .  ,  e  1  =  0 
L  jj,L  ) 


KTTl  „  -  -KTTli 

V<A  T_.  ,e  ,e  1=0 

V  ‘  L  oj  ,  L  L  oj  ,L  / 


(28) 


(29) 


where  j  =1,2,3.  Equations  (28,29)  may  be  written  in  matrix  notation  as 

K(K)q  =  0  (30) 

where  q  is  a  column  matrix  whose  elements  are  A  and  B  (L= 1,2,3), 

“»>  L  L 

while  K  is  a  6x6  square  matrix  and  is  a  function  of  k.  For  a  nontrivial 
solution  of  q  we  must  have 

l|K(<3|l=0  (31) 


which  provides  the  roots  k. 

KTri  “iCTTi 

Notice  that  when  <  =  4,  e  =  -e  =  i  and  the  left-hand  side  of 
Eq.  (28)  is  identical  to  that  of  Eq.  (29)  except  the  minus  signs.  Hence  the 
first  three  rows  of  K  are  identical  to  the  last  three  rows  except  the 
signs  and  the  rank  of  K  is  no  more  than  3.  This  means  that  k  =  4  is  a 
root  of  multiplicity  of  at  least  3  and  that  there  are  at  least  3  independent 
solutions  for  q  in  Eq.  (30)  associated  with  k  =  ^5. 

To  find  q,  i.e.,  A  and  B  for  <  =  h>  let 

L  L 


L 


3j  ,  L  3s ,  L 


J 


(52) 
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where  6.  is  the  Kronecker  delta.  Thus  regarding  x  .  ,  as  a  3x3  matrix, 
J  S  3  J  ,  L 

xt.  is  the  adjoint  matrix  and  J  is  the  determinant.  Assuming  that 
,  L 

J  ^  0,  the  solution  given  by 


Ay  =  —j=-  J‘^k  T*  ,  , 

L  2v2  s  3s, L 

satisfies  Eqs.  (28,29)  when  <  =  %. 
real  constants.  With  Eq.  (33),  Eq. 
that  <  =  h, 


"  2/2 


j-^k  T* 


s  3s,  L 


(33) 


k^  (s  =  1,2,3)  in  Eq.  (33)  are  arbitrary 


(18)  can  now  be  written  as,  noticing 


(r,4)) 


1 

(2r)^^ 


L 


J“‘x* 


T 

3s,L  ij,L^L 


(34) 


Equation  (34)  provides  three  independent  solutions  for  the  stress  dis¬ 
tribution  near  the  crack  tip.  Notice  that  at  <}>  =  0,  Eq.  (34)  for  i  =  3 
has  the  expression,  after  using  Eqs.  (19)  and  (32), 


a  (r,0)  =  ^^^k.  ,  j  =  1,2,3  (35) 

3J  (2r)'^  2 

Hence  k^  are  the  stress  intensity  factors.  It  is  shown  in  the  Appendix 
that  the  solution  associated  with  k^  can  be  identified  with  the  solution 
near  the  tip  of  a  crack  length  2a  due  to  a  uniform  symmetric  tensile  stress 
0^2  at  infinity.  Similarly,  the  solutions  associated  with  k^  and  k^ 
can  be  identified  with  the  solutions  near  the  tip  of  a  crack  length  2a  due 
to  uniform  antisymmetric  shear  stresses  and  respectively,  at 

infinity. 

Equation  (34)  is  not  valid  when  J  =  0.  However,  J  =  0  implies  that 
J  =  0  and  the  rank  of  K  is  no  more  than  two.  It  follows  that  k  =  4  is 
a  root  of  multiplicity  of  at  least  four  when  J  =  0.  This  is  not  possible 
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because  for  the  special  case  of  a  crack  in  an  isotropic  material  it  is  known 
that  K  =  is  a  root  of  multiplicity  three. 


4.  SINGULARITIES  AT  THE  TIP  OF  A  CRACK  NORMAL  TO  AN  INTERFACE 


We  will  now  proceed  to  determine  the  singularities  at  the  tip  of  a 
crack  which  is  normal  to  the  interface  between  two  anisotropic  materials  as 
shown  in  Fig.  2.  Equations  (17)  and  (18)  for  the  displacements  and  the 
stresses  apply  to  materials  1  and  2.  To  distinguish  notations  in  material  2 
from  those  in  material  1,  we  will  add  a  prime  to  all  notations  referring  to 
material  2  except  <  which  is  the  same  in  both  materials  and  r  and  (p 
which  have  unambiguous  definitions.  Since  material  1  is  divided  into  two 
parts  by  a  crack,  the  displacements  and  the  stresses  in  both  parts  need  not 
be  given  by  the  same  expressions.  We  will  therefore  use  a  superscript  (+) 
to  denote  quantities  referred  to  material  1  in  the  x^  >  0  region  and  a 
superscript  (-)  to  denote  quantities  referred  to  material  1  in  the  <  0 
region. 

The  stress-free  boundary  conditions  at  the  crack  surface  are 


+ 

0,  . 

=  0 

at 

([)  =  TT 

=  0 

at 

())  =  -IT 

and  the  continuity  conditions  at  the  interface  are 


u'!’  -  u!  =0 
1  1 

a'^  -a'  =  0 

2j  2j 


u;  -  u.  =0 
1  1 


a'  -  0^.  =0 


at  <J)  =  tt/2 


2j  2j 


(56b) 


at  <j)  =  -tt/2 


(56c) 
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where  i,j  =  1,2,3.  Use  of  the  expressions  for  and  from  Eqs.  (17) 

and  (18)  in  Eqs.  (5ba) ,  (56b),  (56c)  results  in  18  equations  for  the  eighteen 
coefficients  B|^,  A^^  and  B^  (L  =  1,2,3)  which  can  again  be 

written  in  the  form 


K(<)q  =  Q 


(37) 


where  K  is  now  an  18x18  square  matrix  and  the  elements  of  q  are  A  ,  B,  , 

Jj  jj 


A^,  B^,  Aj^  and  Bj^  (L  =  1,2,3).  For  a  nontrivial  solution  of  q, 


we  must 


have 


K(<)  II  =  0 


(38) 


It  should  be  pointed  out  that  for  real  <  one  could  use  the  real 
expressions  for  u^  and  in  Eqs.  (21)  and  (22).  We  again  obtain  18 

equations  in  the  form  of  Eq.  (37)  in  which  the  elements  of  q  are  a*, 

Li 

a^,  a^,  a^,  a^  and  a^^.  The  elements  of  the  18x18  matrix  K  are  real. 

The  roots  <  of  Eq.  (38),  whether  real  or  complex,  can  be  found 
numerically  when  the  elasticity  constants  c. .  and  c! .  of  both  materials 
are  known.  In  the  complex  plane  k,  one  can  find  the  number  of  roots  < 
in  a  given  region  by  using  the  following  theorem  [22]:  If  a  function  >J^(<)  , 
whose  only  singularities  inside  a  closed  contour  C  of  the  K-plane  are 
poles  and  whose  value  is  not  zero  at  any  point  on  the  contour,  then 


A/ 


d:i;/dK 


dK  =  N  -  P 


(39) 


where  N  is  the  number  of  zeros  and  P  is  the  number  of  poles  inside  C. 
In  our  problem  ip  =  ||k||  and  Eq.  (39)  can  be  written  as 


^  /  ‘'1'" 


K  =  N 


(40) 


because 


11k  11  is  bounded  and  hence  there  is  no  pole.  In  view  of  the  fact 


I  that 

!ln  4)  =  Jin  14^  1  +  i  arg  4^  (41) 

Equation  (40)  reduces  to 

'  N  =  ^  (arg  lljCll)  (42) 

where  denotes  the  change  in  the  value  of  f  in  going  around  the  con- 

I  tour  C.  Equation  (42)  implies  that  if  1|k1|  has  N  roots  inside  the 

contour  C  in  the  <-plane,  the  locus  of  HjCH  when  plotted  on  the  complex 
lljCll  plane  encircles  the  origin  N  times  as  <  goes  around  the  contour  C. 

I  This  is  a  useful  result  which  can  be  employed  to  determine  the  number  of 

roots  in  a  given  region. 


5.  NUMERICAL  EXAMPLES 


For  numerical  examples  we  will  assume  that  each  layer  in  the  composite 
shown  in  Fig.  1  is  an  orthotropic  material  with  respect  to  the  (x^^,  x^,  x^) 
axes  where  x^  =  x^  and  the  x^-axis  is  the  fiber  direction  which  makes  an 
angle  0  with  the  x^-axis.  We  further  assume  that  all  layers  are  identical 
orthotropic  materials  although  the  ply-angle  9  may  vary  from  layer  to 
layer.  Referring  to  the  (x^^,  x^,  x^)  axes,  the  following  engineering  con¬ 
stants  for  two  different  orthotropic  materials  are  taken  from  [23]  and  [24], 
respectively. 

Composite  W 

(Typical  high  modulus  graphite/epoxy,  [2.3]) 


E  =  E^  =  2.1  x  10  psi 
12  ^ 


=  20  X  10  psi 


(43) 


Gi2  =  0^2  =  =  0.85  X  10  psi 


''21  '  ^31  =  '’32  '  “•2' 


I 


I 


i 


I 


I 


I 


I 


I 


Composite  T 

CT300/5208  graphite/epoxy,  [24]) 


E  =  E  =  1.54  X 10  psi 


E_  =  (20  X 10  psi)' 


Gi2  =  G23  =  G,,  =  0.81x10  psi 


31 


'’21  ’  ''31  =  ''32  = 


(44) 


In  Eqs.  (43)  and  (44),  E^  are  the  Young's  moduli,  the  shear  moduli 

and  are  the  Poisson's  ratios  [25].  They  are  related  to  [25,26], 

and  hence  s. .  can  be  computed.  Since  E. .,  G. .  and  v. .  given  in 
Eqs.  (43)  and  (44)  are  referred  to  the  (x  ,  X-,  x  )  axes,  s..  so  obtained 

X  4w  1 J 

will  be  denoted  by  s^^ .  Its  inverse  c^^  is  obtained  by  using  the  equa¬ 
tions  derived  in  [25] .  Finally,  we  obtain  and  c^^  which  are 

referred  to  the  (x,  ,x-,x_)  axes  from  s.  .  and  c.  .  (see  [2,12])  when  the 
ply-angle  0  in  material  1  is  given.  For  a  different  ply-angle  0'  in 

material  2,  we  obtain  s!.  and  cf .  in  the  same  way. 

13  13 

With  the  elasticity  constants  c.  .  and  c.'.  associated  with  the  ply- 

13  13 

angles  0  and  0'  in  materials  1  and  2  so  determined,  the  formulation  in 
Section  4  provides  the  matrix  K  and  the  roots  of  the  determinant  ||k|1 
furnishes  the  desired  roots.  In  Table  1  we  list  the  roots  K  for  the 
Composite  W  given  by  Eq.  (43)  for  combinations  of  (0/0')  angles  in  which 
both  0  and  0'  assume  the  values  between  -90°  and  90°  with  an  increment 
of  15°  degrees.  Since  the  values  of  k  for  (-a/b)  and  (a/-b)  composites 
are  identical,  it  is  sufficient  to  consider  0  >  0  only.  In  Table  2  we 


"'The  value  of  E3  in  [24]  is  22x10^  psi.  The  calculations  presented  in 
Table  2  and  Figs.  4  and  6  were  inadvertently  based  on  E3 = 20  x 10^  psi.  How 
ever,  a  few  sample  calculations  with  the  correct  E3  show  that  the  errors 
on  <  are  no  more  than  1%. 
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list  the  roots  <  for  composite  T.  In  both  tables,  k  =  ^5  is  a  triple 
root  when  9=6'  and  when  9  =  -0’  =  90®.  This  is  expected  since  when 
0=9'  and  9  =  -0'  =  90®  material  1  and  material  2  have  the  same  fiber 
orientation  and  hence  there  is  no  real  interface  between  material  1  and 
material  2.  Thus  the  crack  is  in  a  homogeneous  material  and  the  analysis 
in  Section  3  shows  that  k  =  %  is  a  triple  root.  For  all  other  combina¬ 
tions  of  (9/9')  there  are  three  non-equal  roots  of  positive  k. 

Notice  that  the  values  of  k  are  identical  for  (0/90)  and  (9/-90) . 
Likewise,  k's  are  identical  for  (0/9')  and  (0/-9')  as  well  as  for 
(90/9')  and  (90/-9').  Some  of  the  roots  are  not  easy  to  locate.  For 
instance,  two  of  the  three  roots  for  (0/90)  in  Composite  T  differ  by  only 
10  At  the  first  calculation  we  only  located  the  root  k  =  0.697726  and 
completely  missed  the  other  two  roots.  Great  effort  and  care  were 
exercised  in  locating  the  two  roots  which  are  very  close  to  each  other. 

We  also  used  Eq.  (42)  to  verify  numerically  that  there  are  indeed  three 
positive  roots  for  k. 

The  results  presented  in  Tables  1  and  2  are  shown  in  graphs  in 
Figs.  3-6.  In  Fig.  3  contour  lines  for  the  largest  c,  denoted  by  k^, 
are  shown  for  Composite  W  and  the  corresponding  contour  lines  for 
Composite  T  is  shown  in  Fig.  4.  In  Figs.  5  and  6  the  contour  lines  for 
(Kj^  -  Kj)  where  is  the  smallest  of  the  three  roots  are  shown  for 

Composite  W  and  Composite  T,  respectively.  In  view  of  symmetry,  only 
half  of  the  contour  lines  are  shown  in  Figs.  3-6.  In  both  composites 
we  see  that  the  largest  occurs  at  (0/-90)  which  is  identical  to 

(0/90).  At  (0/90),  =  0.676635  for  Composite  W  and  =  0.697726 

for  Composite  T.  The  largest  difference  between  and  occurs 

approximately  at(30/-40)  for  both  Composite  W  and  Composite  T. 
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6.  UNIFORM  EXTENSION 

We  will  consider  in  this  section  the  problem  in  which  the  composite 
shown  in  Fig.  2  is  subject  to  a  uniform  extensional  strain  in  the 

x^-direction.  For  this  purpose,  we  add  the  term  to  the  right-hand 

sides  of  Eqs.  (5),  (13),  (17)  and  (21)  for  the  displacements  u^  and  the 
term  C7),  (14),  (18)  and  (22)  for  the  stresses  [1,27]. 

With  these  modifications,  application  of  the  boundary  and  interface  condi¬ 
tions,  Eqs.  (36),  results  in  the  following  system  of  equations 

r‘'^K(<)q  =  £jb  (45) 

where  the  elements  of  q  are  A^,  B*,  A|^,  B|^,  A~  and  B^  while  b  con¬ 
tains  Since  the  right-hand  side  of  Eq.  (45)  is  independent  of  r, 

we  let  K  =  0  and  obtain 

K(0)q  *  £^b  (46) 

Equation  (46)  provides  q  if  K(0)  is  not  singular.  If  K(0)  is  singular, 
a  solution  may  still  exist  if  b  is  orthogonal  to  the  left  eigenvectors  of 
K(0).  However,  when  k  =  0  Eq.  (18)  indicates  that  are  constants  and 

the  stresses  are  uniform.  Therefore,  instead  of  solving  Eq.  (46),  we  will 
use  an  alternate  approach  employed  in  [1]  to  find  a  uniform  stress  solution 
which  satisfies  the  prescribed  uniform  extensional  strain 

To  this  end,  we  first  solve  from  the  first  of  Eq.  (25) ; 


(j  1) 


We  then  eliminate  a,  from  the  rest  of  Eq.  (25)  to  obtain 


e.  =  R. .a.  +  R, e,  , 
^  11' 


(i  j  1) 
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while  the  expressions  for  e,  and  in  material  2  are 

£_  =  +  R' 

j  o2  2  j3  o  35  5  jo  6  3  1 

(54] 

=  R'  a  + R'  a'  + R'  a'  + R'  a  + R'e 
5  52  2  5j  3  55  5  56  6  51 

Again,  the  superscripts  +,  -  and  a  prime  for  £„  and  are  omitted  in 

view  of  Eq.  (51]. 

For  the  composite  considered  in  Eq.  C43)  and  (44],  the  material 
property  in  each  layer  is  symmetric  with  respect  to  the  *  0  plane  and 
hence  [24,25], 


s,  ,  =  s,,  =  s«.  =  S-,  =  s,.  =  s_,  =  s_.  =  s_,  =  0 

14  16  24  26  34  36  54  56 


Tliis  implies  that  R_^  =  R'  =  R_,  =  R'  =  0  and  elimination  of  £_  and 

36  36  56  56  3 

£g  from  Eqs.  (53]  and  (54)  reduces  to 

(56) 

Equations  (56]  provide  a  one-parameter  family  of  solutions  for  a^,  and 
is  arbitrary  while  a^,  a”  and  are  determined  from  Eq.  (47). 

This  completes  the  solution  for  non- zero  elements  of  the  stresses  shown  in 
Eq.  (52).  It  is  not  difficult  to  see  from  Eqs.  [52]  and  (47)  that  ^7^  =  . 


f''32-'‘32>°2*'>33“3*’>35“5  ' 


Equations 

(56) 

have  no 

solutions 

if 

*^32 

"35  , 

"3  ~  "3 

^^52  - 

"52 

"53 

'  ^5 

"5  -  "5 

(57) 


For  the  composites  considered  here,  Eqs.  (57)  do  not  hold  and  a  uniform 
stress  solution  exists  for  a  given  extensional  strain  Hence,  there  is 

no  additional  stress  singularity  at  the  crack  tip  due  to  the  uniform 
extension  for  the  composites  considered  here.  If  Eqs.  (57)  held,  one 
would  have  a  logarithmic  stress  singularity  at  the  crack  tip  H- 


CONCLUDING  REMARKS 


/  . 

IVe  have  shown  that  the  stress  singularity  at  the  tip  of  a  crack  which 

is  normal  to  and  ends  at  an  interface  between  two  anisotropic  elastic  media 

has  the  form  r  where  <  is  a  root  of  Eq.  (38).  We  show  that  there  are 

in  general  three  distinct  positive  roots  of  <,  say  <  (s=l,2,3).  For  each  < 

s  s 

Eq.  (37)  furnishes  the  associated  q  whose  elements  are  A  ,  B  in 

Eq.  (18).  Since  q  is  unique  up  to  a  multiplicative  constant,  say  c  , 

•<-S  s 

the  solution  given  by  Eq.  (18)  for  <  =  k.^^,  and  can  be  superimposed 

and  written  as 


a.  . 
tj 


•^1  (1) 
c,r  ^  a.  . 

1  ij 


-< 


+  c.,r 


ij 


■'^3  (3) 
+  c  r  a 
3  ij 


(58) 


C  s ) 

where  (s= 1,2,3)  are  known  functions  of  $.  The  analysis  presented 

f  s  1 

here  provides  <  and  a.  .  (s  =  1,2,3).  The  unknown  constants  c.  ,  c- 

S  1.  J  X.  tu 

and  c,  have  to  be  determined  from  the  complete  boundary  conditions  of  the 
whole  composite.  One  could  use  the  solution  obtained  here  to  form  a  special 
finite  element  at  the  crack  tip  and  use  regular  finite  elements  elsewhere  to 
find  a  complete  solution  numerically  for  the  stress  distribution  in  the  com¬ 
posite  with  a  crack. 

It  should  be  pointed  out  that  Eq.  (38)  provides,  besides  the  three 
positive  <  presented  here,  infinitely  many  negative  ic  whose  real  part  is 
negative.  Thus  to  increase  the  accuracy  in  the  special  finite  element  at 
the  crack  tip,  one  may  add  as  many  terms  as  one  pleases  to  the  right-hand 
side  of  Eq.  (58).  These  additional  terms  are  non-singular  at  the  tip  of  the 
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APPENDIX 


Consider  a  crack  of  length  2a  located  on  the  x^-axis  between  x^  =  -a 
and  x^  =  a  in  an  infinite  anisotropic  elastic  material  subjected  to  a 

QO  00 

uniform  tensile  stress  as  well  as  a  uniform  plane  shear  stress  c, 

oo  o2 

OO 

(,Fig.  7)  and  a  uniform  antiplane  shear  stress  (not  shown  in  Fig.  7) 

at  infinity.  The  surface  of  the  crack  is  stress-free.  This  problem  was 
studied  in  [l6,17]  in  which  the  material  property  is  assumed  to  be  symmetric 
with  respect  to  the  x^  =  0  plane  and  in  [l8,19]  for  general  anisotropic 
materials.  We  will  present  an  alternate  expression  for  the  solution  in  the 
entire  plane,  not  just  near  the  crack  tip  as  in  [l9].  For  the  present 
problem,  it  suffices  to  consider  the  associated  problem  in  which  the  crack 
surface  is  subjected  to  the  stresses 


a, .  =  -0, .  on  -a  <  Xt  <  a 

jj  oj  2 


(Al) 


and  the  stresses  vanish  at  infinity.  To  this  end,  we  choose  in  Eqs.  (.13,14) 


-  I  /T 

-v'- 


(A2) 


where  are  arbitrary  complex  constants, 
the  displacements 


We  then  have  the  solution  for 


u.  =  ^  Re 


(A3) 


and  the  stresses 


a.  .  =  y*  Re 


A,  T ...  1  -  -  r 

L  IJ.L  \(:, -.-D  ^  y 


(A4) 


which  vanish  at  infinity.  To  determine  A^,  we  apply  the  boundary  condi¬ 


tions,  Eq.  (Al) : 


(  /  -^2 

-a..  =  •  r  1  2  2^='  ^ 

Jj  L  jj,L  y±i(a2-x|)* 


,  -a  <  X-  <  a 


Equation  (AS)  is  satisfied  if  we  let 


LOO 

I,.  ,A,  =  a,. 
jj.L  L  oj 


1  ^  * 

A.  = 

L  3s  js,L 

where  t*  is  defined  in  Eq.  (32).  We  therefore  obtain  the  solution 
3s ,  L 


a.  .  Re 


a.  .  =  5^  Re  r"^-  -  r  I 

Ij  ^  os.L  1),L  \  (22  .  a2) 


-4 


To  find  the  stress  distribution  near  the  crack  tip,  say  x^  * 


notice  that  as  Zj^  -  a  =  0 


(Zl-  a)  ^-a 


(AlO) 


If  we  let  Z.  -  a  =  tq  and 
L  L 

k  =  /i:a! 

s  os 


(All) 


Equation  (A9)  reduces  to  Eq.  (34).  Thus  the  solution  associated  with  k.  in 

00 

Eq.  (34)  is  due  to  the  uniform  symmetric  tensile  stress  a,„  applied  at 
infinity  while  that  associated  with  k^  and  k2  are  due  to,  respectively, 

30  CO  , 

the  uniform  anti -symmetric  shear  stresses  and  applied  at 

00 

infinity.  It  should  be  pointed  out  tnat  the  solution  associated  with  a^- 

is  not  necessarily  symmetric  with  respect  to  the  X2  axis  even  though  the 

00 

loading  at  the  infinity  is.  Likewise,  the  solutions  associated  with 
and  0*2  are  not  necessarily  anti-symmetric  with  respect  to  the  X2  axis. 
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Fig.  2  A  crack  which  is  normal  to  and  ends  at  the 

interface  between  two  anisoptropic  materials 
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